4! 


*  *'  *4'  •*  ;.,i 


a«sS£;#fe  ©ssiteSf' 

..  ■  .  :.  .»  _  .*  '  <s»ri“  « .:TA. .  <* 


,  •  -  -  ^  *<, 

L.  'ir*!»r-?--~  -» «* 

-  -  *;  ■  -  >  ^  V*  A 


„  u^m 

*(  ^  >■;  * 


"  V  %  :.;  lT5^  ‘  • 

■'  ■■'%.' Jr. »  o.  .^v  <  > 

„./'  ' ",  &"*  U/J  It  '  :>^  ^ 

..  >/.%■■!.' 

'  V.  <,;'-'*-'■■■■■ '  '.■ 

•;  *?  ■  *  -  *  ;  ■;■  ,  „  • 
::■*  "•  •  .  /■..'  ^  ,  V 

<r-s  A  -  o 

<*  >  ; *  V  **  •  ■<  .-■  j 


C  t  E  A  5  {  N  8  H  0  0  S  C 

frugal  sastmnc  anj> 

TgCH^CAf-  ‘  •JFCu'IATIOH  _ 

ijljA  \  Ji  ' 

iif*  mm 


>:y  y-  ‘  ,,.  .'  *  '.. ./  ■  ■  •,  "  .@*~i d&L'  [ 

ft  '*/ u  <3  '*  ,  v.  •'  .  ^  ‘  S. 

Stop  Load  Moving  with  Superseisjnic  Velocity  on  the 
: ■  '  Surface  of  ah  Elastic-Plastic  Half-Space 


•  ky 


■0,  ,  o' 


1  •  .  .  ••  . 

‘  H.’  H*  Bleich  .  and  K>  T .  Matthews 


.  VVv  - 


'  *•  ""V": 


StelrO- 


k  -  v^j 

5-  _•***! 


-  ■  •■•  -  • 


-  Office  of  Naval  Research 
Project  NR,  064“-428>> , 

.  _ -Nohi’*— ~* 

^  I’echnipal  Report  -  .No ,  38 

-  f  dNR-266  (86 )  -CE 

V,-1' ^  J  .  V 

-  Decentber  .1965 

i.  '  1  *  '  .  T 


■  P-Wciluctl<3fl  if!  vhalt  cr  lYi  fl.-i  is  for  any  parfs>#« 

»i  thfc  rfnited  Makui  fotff  ru.Scnt 


;■ 


«9ssfe$sr 


Columbia  (Hnibersitp 

mtijeCitpofJtogorfi 

DEPARTMENT  OF  CIVIL  ENGINEERING 
AND  ENGINEERING  MECHANICS 


Step  Load  Moving  with  Superseismic  Velocity  on  the 
Surface  of  an  Elastic-Plastic  Half-Space 


by 


H,  H 0  Bleich  and  A.  Tc  Matthews 


Office  of  Naval  Research 
Project  NR  064-428 
Contract  Nonr-266(86) 
Technical  Report  No.  38 
CU-4-65  ONR-266(86)-CE 

December  1965 


Reproduction  in  whole  or  in  part  is  permitted  for  am  purpose 
of  the  United  States  Government 


t 


y 


ABSTRACT 


The  plane  strain  problem  of  a  step  load  moving  with 
uniform  superseismic  velocity  V  )>  cp  on  the  surface  of  a 
half-space  is  considered  for  an  elastic-plastic  material 
obeying  the  von  Mises  yield  condition. 

Using  dimensional  analysis  the  governing  quasi-linear 
partial  differential  equations  are  converted  into  ordinary 
nonlinear  differential  equations  which  are  solved  numerically 
using  a  digital  computer.  To  overcome  computing  difficulties 
asymptotic  solutions  are  derived  in  the  vicinity  of  a  singu¬ 
lar  point  of  the  differential  equations . 

Typical  numerical  results  are  presented  for  selected 
values  of  significart  non-dimensional  parameters,  i.e.  of 
the  surface  load  pQ/k,  of  Poisson’s  ratio  v  ,  and  of  the 
velocity  ratio  v/cp. 
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LIST  OF  SYMBOLS* 


a 

1  •  • 

•  4 

Functions  defined  by  Eqs.  ( 3— 39 ) * 

b 

1  •  • 

•  4 

Functions  defined  by  Eqs.  ( 3-1^) -( 3“l6) 
and  ( 3-26) . 

V 

CS' 

c 

Velocity  of  propagation  of  elastic  P-waves, 
S-waves  and  of  inelastic  shock  fronts, 
respectively. 

F 

Plastic  potential,  Eq.  (2-1). 

G 

Shear  modulus. 

J  , 

l 

J' 

2 

Invariants  of  stress. 

k 

Yield  stress  in  shear,  Eq.  (2-1). 

K  = 

2(l+v) 

3(l-2v)  G 

Bulk  modulus. 

L  >  0 

Function  related  to  inelastic  behavior, 

Eq.  (2-36). 

p(x- 

-vt) 

Surface  pressure. 

PE' 

PL' 

Po 

Intensities  of  step  load  surface  pressure. 

R 

Variable  defined  by  Eq.  (2-19). 

s  ; 

i 

s 

2 

Principal  stress  deviators. 

S  ) 

X 

V 

s .  , 
i: 

Stress  deviators  with  respect  to  axes 
x, y,  etc. 

t 

Time. 

v  , 

X 

V 

V  ,  V 

x’  y 

Components  of  particle  velocity  and 
acceleration  in  x  and  y  directions, 
respective ly. 

V 

Velocity  of  surface  presssure. 

*) Other  symbols,  which  are  used  in  one  location  only,  are 
defined  as  they  occur. 
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Values  of  X  at  P-  and  S-fronts,  respectively. 
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Small  quantity  for  purposes  of  asymptotic 
expansion. 

Increments  of  a,v,x,  etc.  at  a  front. 

Small  quantity  for  purposes  of  asymptotic 
expansion. 

Strains,  strain-rates. 

Small  quantity  for  purposes  of  asymptotic 
expansion. 

Angle  defining  direction  of  the  principal 
stress  a  ,  Fig.  4. 

Function  related  to  inelastic  behavior, 

Eq,  (2-8) . 

Poisson*  s  ratio. 

Variable  defined  by  Eq.  (2-13). 

Mass  density  of  medium. 

Stresses,  stress  rates. 

Principal  Stresses. 

Shear  stress. 
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INTRODUCTION. 


The  two  dimensional  steady-state  problem  considering  the 
effect  of  a  pressure  pulse  p(x  -  Vt)  progressing  with  the  velocity 
V  on  the  surface  of  an  elastic  half-space,  Fig.  1,  has  been 
treated  by  Cole  and  Huth  [1]  for  a  line  load.  By  superposition 
their  solution  may  be  used  to  find  the  effect  for  any  other 
distribution  p(x  -  Vt) .  The  equivalent  problem  for  linearly 
viscoelastic  materials  was  treated  by  Sackman  [2],  and  Workman 
and  Bleich  [3],  in  the  superseismic  and  subseismic  ranges, 
respectively.  The  present  paper  considers  the  problem  again, 
but  in  an  elastic-plastic  material  subject  to  the  von  Mises 
yield  condition.  In  this  material  thf  problem  becomes  non¬ 
linear,  so  that  superposition  is  not  permitted  and  each  pressure 
distribution  p(x  -  Vt)  poses  a  separate  problem.  The  present 
paper  treats  only  the  case  of  a  progressing  step  load 

p(x  -  Vt)  =  p  Et(vt  -  x)  . 
o 

Interest  in  steady-state  problems  is  not  only  due  to  the 
fact  that  they  are  natural  stepping  stones  towards  physically 
more  meaningful,  but  more  complex  nonsteady  problems.  In  the 
elastic  case  Miles  [4]  has  considered  the  three  dimensional 
problem  of  loads  with  axially  symmetric  distribution  p(r,t) 
over  an  expanding  circular  area  on  the  surface,  Fig.  2.  He 
has  demonstrated  that  the  plane  problem  [1]  is  the  asymptotic 


solution  for  the  three  dimensional  one  in  the  region  near  the 
wave  front.  This  gives  rise  to  the  expectation  that  the 
\  situation  in  other  materials  may  be  similar.  This  motivation 

for  the  search  for  steady-state  solutions  limits  interest  to 
those  which  do  not  violate  conditions  which  asymptotic  solutions 
of  three  dimensional  problems  must  satisfy. 

The  conditions  to  be  imposed  on  steady-state  solutions 
to  eliminate  any  which  are  not  asymptotic  solutions  of  the 
problem  in  Fig.  2,  or  of  a  similar  one,  can  easily  be  recognized 
in  the  elementary  example  of  a  half-space  of  an  inviscid  com¬ 
pressible  fluid  loaded  by  a  uniform  pressure  pulse  p,  which 
progresses  with  supersonic  velocity,  V  c.  There  is  an  obvious 
solution,  Fig.  3a,  in  which  the  load  produces  a  plane  wave  of 
intensity  p  progressing  with  a  front  inclined  at  the  appro¬ 
priate  angle  f  =  sin-1  ^  .  However,  this  is  not  the  only 
steady-state  solution.  An  alternative  is  a  plane  wave,  the 
front  of  which  is  inclined  at  the  angle  180°-'!//.  Combinations 
of  the  two  solutions  are  also  correct  steady-state  solutions. 

To  find  states  generated  by  the  application  of  a  progressing 
pressure  on  the  surface  only,  it  can  be  reasoned  that  solutions 
which  include  the  wave  front  shown  in  Fig.  3b  can  not  apply 
because  the  medium  ahead  of  the  front  shown  in  Fig.  3a  should 
be  undisturbed  when  the  applied  load  advances  from  the  left 
with  supersonic  velocity.  Further,  the  transient  supersonic 
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solution  being  irrotational  without  a  singularity  in  pressure 
or  velocity  at  the  wave  front,  the  same  must  hold  in  the  steady- 
state.  This  reasoning  leads  to  unique  steady-state  solutions 
for  fluids  and  also  for  elastic  materials  in  supersonic  and 
superseismic  cases,  respectively.  An  equivalent  approach  will 
be  utilized  in  this  paper  . 

In  elastic-plastic  materials  different  sets  of  differen¬ 
tial  equations  apply  in  the  elastic,  or  neutral,  and  in  the 
plastic  regions.  The  fact  that  these  regions  have  moving,  a 
priori  unknown  boundaries,  complicates  the  solution  of  dynamic 
problems  considerably.  In  the  following,  the  basic  equations 
will  be  formulated  and  solved  separately  in  plastic  and  in 
nondissipative  regions.  The  partial  solutions  will  be  matched 
to  obtain  a  complete  one  satisfying  the  prescribed  surface 
conditions  and  additional  ones  obtained  from  the  requirement 
that  the  steady-state  solutions  should  qualify  as  asymptotic 
for  a  transient  problem  of  the  type  shown  in  Fig.  3.  Using 
dimensional  analysis,  similar  to  the  approach  used  in  a  simpler 
case,  [5]  permits  transformation  of  the  original  nonlinear 
partial  .ifferential  equations  in  plastic  regions  into  a  set 

"*) 

In  subsonic  or  subseismic  cases  the  equivalent  approach  is 
not  fully  successful.  For  example,  in  elastic  materials  the 
steady-state  solution  in  the  subseismic  range  xs  unique  for 
many,  but  not  for  all  quantities.  Expressions  for  the  hori¬ 
zontal  stress  and  for  the  velocities  contain  arbitrary  con¬ 
stants  which  can  no'-  be  determined. 
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of  simultaneous  ordinary  ones.  Their  solution  in  the  non- 
dissipative  case  is  elementary.  Those  in  the  plastic  case  are 
nonlinear  and  much  too  complex  for  closed  solution,  but  they 
can  be  solved  numerically  requiring  the  solution  of  transcen¬ 
dental  equations  and  quadratures,  for  both  of  which  digital 
computers  are  employed.  Break  down  of  the  computer  solutions 
in  the  vicinity  of  a  singular  point  of  the  set  of  differential 
equations  made  the  derivation  of  asymptotic  solutions  necessary. 

While  the  solution  of  the  problem  is  obtained  without 
overt  use  of  characteristic  methods,  it  is  actually  dependent 
on  the  hyperbolic  character  of  the  quasi-linear  differential 
equations.  As  a  matter  of  fundamental  interest  it  is  proved 
in  Appendix  B  that  the  steady-state  problem  in  plastic  regions 
is  hyperbolic  for  superseismic  velocities  V  )>  cp  . 


2 .  FORMULATION  OF  THE  BASIC  EQUATIONS. 

Figure  4  indicates  the  half-space  and  a  system  of  Cartesian 

coordinates.  The  x-axis  is  in  the  direction  of  motion  of  the 

step  load,  the  y-  and  z-axes  are  normal  to  the  surface  in  and 

out  of  the  plane  of  the  figure,  respectively.  The  analysis 

considers  the  case  of  plane  strain,  e  =  0,  when  the  velocity 

z 

V  of  the  step  load  is  superseismic ,  i.e.  larger  than  the  largest 
elastic  or  plastic  wave  velocity,  which  is  the  one  of  elastic 
P-waves  in  the  material.  Throughout  the  analysis  it  is  assumed 
that  the  strains  and  velocities  are  small,  so  that  their  higher 
powers  may  be  neglected  in  comparison  to  linear  terms. 

To  describe  the  behavior  of  the  elastic-plastic  material 
the  yield  function  F  is  introduced 


F  =  J» 


2 


-  k2 


(2-1) 


where  j«  j_s  the  invariant 
2 


(2-2 


J»  =  —  s.  .  s . . 
2  2  il  31 


and  the  value  k  )  0  is  the  yield  stress  in  shear. 

The  behavior  of  an  element  of  the  material  is  defined  by 


the  following  statements. 


1.  The  value  of  the  function  F  may  never  be  positive 

F  <  0  (2-3) 

2.  If,  in  an  element  of  the  material  at  a  given  instant, 

F  <  0  (2-4) 

the  rates  of  change  in  stress  and  strain  are  related 
by  the  conventional  elastic  relations. 

3.  However,  if  the  yield  condition 

F  =  0  (2-5) 

is  satisfied,  three  possibilities  exist:  a.,  in  the 
next  instant  of  time  the  material  may  be  in  a  state 
of  plastic  deformation;  b.,  it  may  be  in  a  state  of 
elastic  unloading;  c.,  it  may  be  in  a  neutral  state, 
a.  If  the  material  is  in  a  state  of  plastic  defor¬ 
mation 


F 


0 


(2-6) 
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the  total  strain  rate  will  be  the  sum  of  an  elastic 
and  a  plastic  portion 


e 


ij 


(2-7) 


where  e. .  is  obtained  from 
y-3 

relations,  while,  based  on 
potential. 


the  conventional  elastic 
the  concept  of  a  plastic 


= 


d  F 

<?  a. 
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(2-8) 


J\  which  must  be  positive 


k  >  0 


(2-9) 


is  an  a  priori  unknown  function  of  space  and  time, 
to  be  found  as  part  of  the  solution  of  the  problem. 

b.  In  case  of  elastic  unloading  F  <(  0  holds,  and  the 
elastic  stress-strain  relations  apply. 

c.  In  the  neutral  state  F  vanishes  as  in  case  a.,  but 
neither  energy  dissipation  nor  permanent  defor¬ 
mation  occurs,  and  the  elastic  stress-strain 
relations  apply.  In  the  present  problem  neutral 
regions  of  a  particularly  simple  type  will  be 
encountered  in  which  neither  the  stress  nor  the 
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strain  changes., 


0. 


For  the  purpose  of  this  paper  it  is  convenient  to  combine 
elastic  and  neutral  regions,  which  will  be  called  "nondissipative", 
as  opposed  to  plastic  regions,  where  ~K  />  0  ,  indicating  that 
energy  is  dissipated.  In  the  nondissipative  regions  the  changes 
in  stress  and  strain  are  governed  by  the  elastic  relations, 
while  in  plastic  regions,  Eqs.  (7)  and  (8)  apply.  Formally, 
the  equations  in  nondissipative  regions  can  therefore  be  obtained 
by  substitution  of  1  =  0  into  the  differential  equations  derived 
below  for  the  plastic  regions,  and  by  replacing  the  conditions 
p  =  F  =  0  by  the  inequality  (3). 

Substituting  Eqs.  (7)  and  (8)  and  the  elastic  stress- 
strain  relations  into  the  relation 


1  , 

€ij  “  2  <Vi,j 


V.  .) 


(2-10) 


the  following  constitutive  equations  are  obtained  for  the  case 
of  plane  strain 


1 

.  l-2v 

T 

aVx 

2  G 

sx  +  6(l+v)G 

J 

1 

+  A  S  =  v- — 
x  dx 

1 

•_  .  l-2v 

» 

T 

C)V 

V 

2G 

Sy  6(1+v)g 

J 

1 

+  A  s  =  ~ 
y  oy 

( 2-lla) 
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>  (2-llb) 


2G 

JL_ 

2G 


T  +  7vT  = 


( s  +  s  ) 
x  .  y ' 


_ _ X 


l-2v 


6( l+v)G 


+  *  (sv  + 


x 


S  ) 

y 


0  J 


J  is  the  first  invariant  of  stress,  s  .  s  and  v  ,  v 
x  xJ  y  x  y 

are,  respectively,  the  stress  deviators,  and  the  components 
of  the  particle  velocity  in  the  x  and  y  directions. 
Further,  there  are  two  equations  of  motion 


^sx  1  ^v- 

N 

O  X 


3  Sx 


1  .  dr 
+ 


by  P  d  t 


,  oJ,  >  by 

_y  +  i_ A  +  ^L=n— ^ 
by  3  by  bx  p  bt 


(2-12) 


Equations  (11)  and  (12)  and  the  respective  requirements  on  F 
and  Tv  complete  the  formulation  except  for  initial  and  boundary 
conditions . 

In  the  steady-state  problem  of  a  half-space  subject  to  a 
surface  pressure  p(xjt)  ,  the  latter  and  all  expressions  for 
stresses,  velocities,  etc.,  in  the  solution  must  be  functions 
of 


4  =  x  -  Vt 


(2-13) 


Equations  (11),  (12)  may  therefore  be  reduced  to  a  set  of  partial 

differential  equations  in  the  two  independent  variables  £  and  y 


_V  ^fx  v(l-2v)  _ 

~  2G  hi  "  6(1+v)g  hi 


hJ 


V 

2G 


^ Sy  _  V(l-2v)  ^Ji 
6(  1+v)  G 


+  s  = 
x 


+  k  s  = 


I  Si  ,  , 

2G  ££  +  *T  " 


V 

2G 


o  s 


X 


SI 


+ 


Sv 

Sf 

Sv 


X 


2£ 


Sy 


,  V( l-2v)  °Ji  ,  .  / _ 
+  s  i+v)g  sT  ( 


}  (2-14) 


+ 


s  ) 
y 


=  o 


j 


5sx  ,  1  3J1  ,  St 
34  3  94  oy 


_ y  .  I _ l  .  hi 

Sy  +  3  hy  +  hi 


-  pv 


Ql  -  _  T  T 

~  "  Pv 


Sv 


S| 


SI 


(2-15) 


In  plastic  regions  the  additional  equation  F  =  0  applies, 

so  that  there  are  a  total  of  seven  relations  for  the  seven 

unknown  quant it res  s  ,  s  ,  t/  J  >  v  ,  v  and  k  \  0.  In  non- 

x  y  i  x  y  ' 

dissipative  regions  Eqs .  (14),  (15)  apply,  but  the  function 

7\  vanishes  identically,  k  =  0,  while  F  must  satisfy  either 

« 

F  <(  0,  or  the  two  conditions  F  =  0,  F  <^_  0  simultaneously.  In 
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the  nondissipative  case  there  are  only  six  differential  equa¬ 
tions  and  six  unknown  quantities.  The  complete  solution  of  the 
problem  is  to  be  obtained  from  the  six  differential  equations 
(14) ,  (15)  and  the  applicable  relations  on  7\  and  F,  subject  to 

appropriate  boundary  or  initial  conditions  at  the  surface  and 
at  the  junctions  of  the  as  yet  unknown  regions. 

It  is  demonstrated  in  Appendix  B  that  the  system  of 
differential  equations  in  plastic  regions  is  hyperbolic  for 
the  case  under  consideration,  V  )>  Cp ,  there  being  four  char¬ 
acteristic  directions  which,  due  to  the  nonlinearity,  are 
stress  dependent.  The  problem  in  nondissipative  regions  [1] 
is  also  hyperbolic,  but  without  the  stress  dependency.  With¬ 
out  the  subdivision  of  the  domain  into  regions  of  unknown 
extent,  e.g.  in  the  purely  elastic  superseismic  problem  [1], 
uniqueness  could  be  established  by  available  theorems  on  the 
initial  and  boundary  conditions  required.  In  the  present  case 
of  mixed  regions  of  unknown  extent  no  such  general  theorem  is 
available,  and  the  question  of  uniqueness  of  transient  as  well 
as  steady-state  solutions  can  not  be  answered  in  general.  For 
the  solution  obtained  here  uniqueness  will  be  demonstrated  by 
detailed  examination  of  all  possibilities  for  solutions  which 
satisfy  the  various  conditions  listed  below. 

1.  On  the  surface,  y  =  0,  a  step  pressure  p  =  pQH(vt  -  x) 


T  '  - 

f 
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1 


•l 


normal  to  the  surface  is  applied,  so  that,  with 
reference  to  Eq.  (13) 


f~  PQ  for  i  <  0 

! 

t 

i  0  for  £,  )>  0 


(2-16) 


while 


T  =  0 


(2-17) 


2.  It  is  known  from  a  general  study  of  elastic-plastic 
wave  propagation  [6]  that  the  largest  characteristic 
velocity  possible  is  cp(  the  velocity  of  elastic 
P-waves.  All  stresses  and  velocities  must  therefore 
vanish  outside  the  wedge  formed  by  the  negative  £-axis 
and  the  P-front  which  is  inclined  at  the  angle 


0 


P 


=  7T 


sin 


-l 


I  J  2G(  1-v)  ~| 
V  1/  p(l-2v) 


(2-18) 


with  the  £-axis,  Fig.  5. 

3.  While  the  character  of  the  solution  of  the  transient 
elastic-plastic  problem  near  £  =  y  =  0  is  not  known, 
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in  the  corresponding  elastic  problem  no  singularity 
occurs  in  the  region  of  nonvanishing  solutions  for 
stresses  and  velocities.  Subject  to  later  confir¬ 
mation  that  the  steady-state  problem  permits  solutions 
without  singularities  such  solutions  are  prescribed. 

It  is  subsequently  seen  that  they  exist,  and  are 
unique . 


To  apply  dimensional  considerations  new  variables 


R  =  -yjz 2  + 


(2-19) 


0  =  COt  1  (^) 


(2-20) 


are  introduced.  The  differential  equations  and  the  various 
additional  conditions  which  determine  the  solution  contain 
only  the  dimensional  quantities:  p,  V,  G,  k,  p^  and  the 
coordinate  R.  However,  three  of  these  six  quantities  can  be 
expressed  by  three  others  and  by  three  nondimens ional  para¬ 
meters.  A  suitable  independent  set  of  dimensional  quantities, 
p  ,  V  and  R,  is  used,  while  the  remaining  quantities  enter 

0  p 

pV2  G 

the  problem  only  in  the  nondimensional  combinations  >  £ 

and  .  The  stresses,  velocities  and  the  quantity  A  in 

terms  of  functions  of  nondimensional  variables  become 
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(2-21) 


S.,  T )  J  =  p  f(0, 
1  10 


vt  =  Vf(0, 


(2-22) 


~k 


(2-23) 


The  first  two  relations  satisfy  requirement  3.  that 
stresses  and  velocities  are  not  singular  at  R  =  0  .  The 
quantity  \  has  a  singularity  for  R  =  0  ,  hut  this  is  not 
objectionable  because  only  its  sign,  not  its  value  is  of 
physical  significance.*) 

Using  the  new  variables  R.,  0,  and  noting  that  the 
expressions  (21),  (22)  are  not  functions  of  R,  one  finds  for 

the  derivatives  in  Eqs.  (14),  (15) 


c5  _  sin  0  d 

=  R  30 


(2-24) 


3  _  cos  0  d 

3y  R  30 


In  this  fashion  the  partial  differential  equations  (14),  (15) 

in  £,  and  y  become  ordinary  ones  in  the  variable  0.  Because 

of  the  manner  of  solution  to  be  employed  the  unknowns  s  .  s 

x  y 

*)  in  addition,  as  7\  does  not  occur  in  the  elastic  case  used 
as  guide  one  can  conjecture  that  it  may  show  the  same  behavior 
in  a  transient  problem. 
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and  t  are  replaced  by  three  new  dependent  variables  s  ,  s 

1 

and  Q, 


s  =  s  cos2  0  +  s  sin2  0  (2-26) 

^  1  2 

s  =  s  cos2  0  +  s  sin2  0  (2-27) 

y  2  i 

t  =  (s  -  s  )  sin  0  cos  0  (2-28) 

12  ' 

s  and  s  are  the  two  principal  stress  deviators  while  0 
12 

is  the  angle  between  the  direction  of  s  and  the  horizontal, 
Fig.  4.  Introducing  further  the  angle  7  between  the  direction 
of  s^  and  the  position  vector,  Fig.  4, 


7  -  <t>  -  e 


(2-29) 


the  six  differential  equations,  (14) ,  (15)  become  finally 


s'  cos2  0  +  s*  sin2  0  -  (s  -  s  )  0'  sin  20  + 


+ 


2 

l-2v 

3(l+v) 


J’  +  L(s  cos2  0  +  s  sin2  0)  = 
11  2 


2G  . 
—  v1 
V  x 


(2-30) 


s’  sin2  0  +  s'  cos2  0  +  (s  -  s  )  0*  sin  20  + 

1  2  12 

+  3(1+77  J1  +  L(si  s':n2  e  +  s2  cos20) 


2G  ,  .  , 

—  vj,  cot  0 

(2-31) 
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(2-32) 


W  Ji  ■  f  cot  *  - 


( s*  -  s»  )  sin  20  +  2(  s 
12  1 


-  s  )  0f  cos  20  +  L(s  -  s  )  sin  20  = 


4r  (v*  cot  0  -  v* ) 
V  '  x  v  y 


(2-33) 


s*  cos  0  sin  7  +  s'  sin  0  cos  7  -  (s  -  s  )  0*  cos  (7-0)  + 
1  2  12 

1  . 


+  sin  0  J’  =  -  pv  v1  sin  0 
3  1  K  x 

(2-34) 


s’  sin  0  sin  7  -  s’  cos  0  cos  7  +  (s  -  s  )  0’  sin  (7-0)  + 

1  2  12 

-  -i-  cos  0  J1  =  -  pV  v1  sin  0 

o  1  y 

(2-35) 

Primes  indicate  differentiation  with  respect  to  0,  and  L  is 
related  to  ~h, 


=  2GR 
V  sin  0 


(2-36) 


The  function  L  is  subject  to  the  same  conditions  as  i.e., 
L  )>  0  in  plastic  regions,  L  =  0  elsewhere. 

The  expression  for  the  plastic  potential  in  these  vari¬ 
ables  is 


T?  -  —2 

1}  — 


H  5  S 
1  2 


_  2 
a 

2 


1.2 
—  A. 


/  /"\  <~r  \ 

) 
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In  plastic  regions  Eq.  (6)  requires  F  =  0,  giving  the  additional 
differential  equation 


(2s  +  s  )  s'  +  (s  +  2s  )  s'  =  0  (2-38) 

121122 

The  unknowns  and  v^,  can  be  eliminated  from  Eqs.  (30)- 
(38)  without  differentiation.  Using  the  symbol 

X  =  x(0)  =  sin2  0  (2-39) 

the  operations 

Eq.  (45a)  =  Eq.  (30)  +  Eq.  (31)  -  Eq.  (32)  (2-40) 

Eq.  (45b)  =  -  X  Eq.  (32)  +  sin  0  Eq.  (34)  -  cos  0  Eq.  (35) 

(2-41) 

Eq.  (45c)  =  X  (sin  20  [Eq.  (31)  -  Eq.  (30)]  +  cos  20  Eq.  (33)}  + 

+  cos  (y  -  0)  Eq.  (34)  -  sin  (y  -  0)  Eq.  (35) 

(2-42) 

Eq.  (45d)  =  X  (cos  20  [Eq.  (31)  -  Eq.  (30)]  -  sin  20  Eq.  (33)}  + 

+  sin  (y  -  0)  Eq.  (34)  +  cos  (y  -  0)  Eq.  (35) 

(2-43) 

Eq.  (45e)  =  Eq.  (38)  (2-44) 

lead  to  the  following  set  of  five  differential  equations  for 
plastic  regions: 


1  +  V 


SOLUTIONS  FOR  INDIVIDUAL  REGIONS . 


3  . 


As  a  first  step  towards  the  construction  of  overall 
solutions,  expressions  for  individual  regions  must  be  derived. 
The  latter  will  be  combined  in  Section  4  to  find  the  solution 
for  the  entire  domain. 


a)  Nondissipative  Regions. 

Equations  (2-46)  are  linear  and  homogeneous  so  that  the 

derivatives  of  the  stresses  s'  ,  s»  ,  J«  and  the  value  (s  -s  )0> 

a  2  l  12 

vanish,  unless  the  coefficient  matrix  in  Eqs .  (2-46)  is 
singular,  requiring 


X(1  -2X)  [  1  +  (1  -  2X)(1  -  2 v )]  =  0  (3-1) 


Equation  (1)  has  two  significant  roots, 


1- v 
1  -  2v 


(3-2) 


and 


(3-3) 


i 

» 
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■53*-  -*  7 


and 


a  spurious  one,  X  =  0  .  Substitution  of  the  two  roots 

Xp  and  Xg  into  (2-39)  furnishes  the  two  locations 


0p  =  7r  -  sin 


-i  fP 

V 


(3-4) 


=7 r  -  sin 


-i  fs 

V 


(3-5) 


where  cp,  Cg  are  the  velocities  of  P-  and  S-waves,  respectively. 

In  all  locations  0  ^  0  or  0  the  values  s',  s*  ,  J*  (s  -  s  )6M 

F  b  12  112 

vanish,  so  that  in  nondissipative  regions  the  stresses  must 
remain  constant  except  at  the  locations  0p  and  0  The  latter 
being  the  potential  locations  of  elastic  P-  and  S-shock  fronts, 
respectively,  it  is  known  that  discontinuities  in  stresses  and 
velocities  may  occur  at  these  locations  and  may,  therefore,  be 
part  of  the  complete  solutions  to  be  constructed.  The  follow¬ 
ing  pertinent  details  will  be  required  subsequently. 

(1)  The  P-front. 

Designating  the  discontinuous  changes  in  the  various 
quantities  at  the  front  by  the  symbol  A,  the  disconti¬ 
nuities  in  the  stresses  =  gz  (normal  and  tan¬ 

gential  to  the  front,  respectively)  and  in  the  component 

_ 

It  is  due  to  multiplication  by  sin  0  in  changing  variables 
and  eliminating  the  velocities. 
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of  the  velocity  (normal  to  the  front)  are  proportional 


to  A oN. 


(3-6) 


No  other  discontinuities  can  occur  in  this  location. 

The  changes  Aa^  and  AcT  are  of  course  limited  by 
the  yield  condition  F  <(  0  which  must  be  satisfied  on 
either  side  of  the  front. 

In  the  actual  solution  a  P-front  will  be  encountered 
only  for  the  special  case  where  the  region  ahead  of  the 
front  is  stressless  and  at  rest.  The  normal  to  the  front 
is  then  a  principal  direction  for  the  stresses  behind  the 

7T 

front,  so  that  y  =  0  or  75-.  Selecting  ,  corre¬ 

sponds  to 


7  = 


7T 

2 


(3-7) 


The  value  of  the  other  quantities  of  interest  behind  the 
shock  front  are 

2(l-2v)a1 

p  =  3  ’  3i  =  3(l-v)  ’  s 

subject  to  the  limitation 


(l-2v)ai  (l+v)a1 

3 ( 1— v )  '  Jx  "  1-v 

(3-8) 
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< 


(3-9) 


a 

1 


V~3(l  -  v)k 
1  -  2v 


imposed  by  the  yield  condition. 

( 2 )  The  S-front . 

At  an  S-front  discontinuities  occur  only  in  the  shear 
stress  =  tT  =  T  an<^  t^ie  tan<3ent;‘-a-*-  velocity  vT. 

The  change  in  velocity  is  proportional  to  At. 


(3-10) 


In  addition,  the  yield  condition  F  <(  0  must  again  be 

satisfied  ahead  of  and  behind  the  front. 

The  relations  between  the  state  of  stress  on  either 

side  of  an  s-front  in  terms  of  At  and  of  the  variables 

s  .  s  .6  and  7  can  be  obtained  in  such  routine  manner 
12 

that  only  one  result  actually  used  in  Section  4  is  pre¬ 
sented  here. 

It  is  possible  for  an  S-front  to  occur  between  two 

neutral  regions,  i.e.  regions  of  constant  stress  for  both 

of  which  the  yield  condition,  F  =  0,  is  satisfied.  in 

this  special  case  the  quantities  (3  ,  J  ,  s  and  s  have 

11  2 

no  discontinuity  at  the  front,  only  the  direction  of  the 
principal  stress  changes.  The  values  of  the  angles  7,  7 
ahead  of  and  behind  the  front,  respectively,  are  comple- 
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mentary 


7 


7T  -  7 


(3-11) 


as  shown  in  Fig.  6  . 

b)  Plastic  Regions. 

In  such  regions  Eqs .  (2-45)  apply.  They  are  linear  and 

homogeneous  in  the  values  s'  ,  s'  ,  etc.,  and  may  be  satisfied 

1  2 

by 

s'  =  s'  =  J'  =  (s  -  s  )0'  =  L  =  0  (3-12) 

1  2  11  2 

However,  L  =  0  implies  A  =  0  which  violates  Eq.  (2-9). 

It  follows  that  in  plastic  regions  the  determinant  of  Eqs. 

(2-45)  must  vanish,  giving  the  "determinantal  equation" 

(s  +  s  )2(b2  +  b  b  )  =0  (3-13) 

1  2  2  1  3 


where 


b  =  2  [1  +  ( l-2v) ( 1~2X) ] 


(3-14) 


b£  =  p  cos  2y  +  (l-2X)(l-2v)  (3-15) 
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(3-16) 


b  =  (l+v)(l-2X)  -  02X 

3 


and 


0 


s  +  s 
1  2 


(3-17) 


If  Eqs.  ( 14)  to  (17)  are  substituted  into  Eq.  (13)  it  is  seen 

that  the  latter  is  a  homogeneous  quadratic  expression  in  s^ 

Due  to  the  vanishing  of  the  determinant  only  four  of  the 

five  Eqs.  (2-45)  are  independent.  By  definition  L  must  not 

vanish,  so  that  s'  ,  s'  ,  J1  and  9 1  can  always  be  expressed 

12  1 

in  terms  of  L, 


s'  = 
1 


(3-e)  b4  (s,  +  s2)  L 


3b. 


(3-18) 


s '  = 
2 


(3+0)  b4  (st  +  ss)  L 


3b. 


(3-19) 


6  1  •■= 


+  sin  2y  bg  L 
0 ( 1-2X)  b  " 


J' 

1 


( 1+v 
2v 


bE  3  ®  b4 


(S1  +  S2)  L 


(3-20) 


(3-21) 


Velocities  and  accelerations  may  be  obtained  from  the¬ 
re  lat ions 
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'7 


-  V  sin  0(s  +  s  )L 

vi  =  - '2G[l^ai:J”b - tb/  sin (2y-*)-  2b3sin  0]  (3-22) 


-  V  sin  0(  s  +  s  ) L 

Vy  =  - 2G(1-2X)  ba  fb2P  cos ( ^Y~0)+  2b3cos  0]  (3-43) 


vx  =  ^  sin  0  (3-24) 

v  =  sin  0  v*  (3-25) 

y  k  y 

where 


b  =  ( 1+  v )  cos  2y  +  (3X(l-2v) 

4 


(3-26) 


Since  Eq.  (13)  must  remain  valid  throughout  a  plastic 
region,  it  may  be  differentiated  with  respect  to  0  .  This 
leads  to  an  expression  which  contains  the  first  derivatives 
of  the  stresses  linearly,  so  that  substitution  of  Eqs .  (18)- 
(20)  furnishes  a  linear  equation  for  L.  Its  solution  gives 
L  as  a  function  of  (3  ,  7  and  of  the  position  angle  05 


3bs(l-2x)Jx  sin  20[4(l-2vXb2+b3)+b1O2+2+2v)]+4b28sin  2y  sin2  0^ 
4  sin2  0  \  (3+p2) (l-2X)b  (b  cos  2y  -  b  $x)  + 3b  b  sin2  2y  I 

V,  4  2  1  2  3  j 


(3-27) 
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The  values  of  the  derivatives  s'  .s'  ,  J'  and  0'  can 

12  1 

be  obtained  by  substitution  of  Eq.  (27)  into  Eqs .  (18)— (21) - 

In  principle  Eqs.  (18)- (27)  permit  the  numerical  deter¬ 
mination  of  the  values  of  stresses  and  velocities  in  the 
interior  of  a  plastir  region  by  quadratures  if  the  values  on 
one  boundary  of  this  region  are  known.  The  starting  values 
must  inherently  satisfy  the  yield  condition,  F  =  0  ,  and  the 
determinantal  equation  (13) .  Further,  the  condition 


L  >  0 


(3-28) 


must  be  satisfied,  to  assure  that  the  result  applies. 


c)  Discontinuities  (Shock  Fronts) . 

It  is  known  that  m  transient  problems  one,  but  just  one 
type  of  plastic  shock  front  can  propagate  in  the  elastic-plastic 
material  considered  here  [6]  .  However,  such  a  front  can  exist 

only  in  locations  where  the  normal  to  the  front  lies  in  the 
direction  of  one  of  the  principal  stresses,  while  the  other 
two  are  equal,  and  where  the  yield  condition  is  satisfied. 

The  velocity  of  propagation  of  the  front  is 


c 


(3-29) 


-26- 


where  K  = 


2( 1+v) 

3  l-2v 


G  is  the  bulk  modulus.  The  discontinuity 


is  restricted  to  the  particle  velocity  normal  to  the  front 

and  to  the  first  invariant  J1  .  The  change  AJ1  must  have 
the  same  sign  as  J 


AJ. 

—  >  0 
1 


(3-30) 


The  other  conditions  stated  define  the  values  of  7  and  3 


7  =  %  ,  3  =  3 


(3-31) 


A  discontinuity  traveling  in  real  space  with  velocity  c, 
Eq.  (29),  can  occur  in  the  steady  state  problem  only  in  the 
location 


0  =  7r  -  sin"1  (£) 


(3-32) 


The  corresponding  value  of  X  is 


(3-33) 


The  denominator  in  Eq.  (2/)  vanishes,  as  expected,  for  these 
values  of  y,  3  and  X. 

The  possibility  of  the  occurrence  of  this  discontinuity 
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-<!  V 


(shock)  must  be  considered  when  constructing  the  complete 
solutions  in  Section  4.  It  was  actually  found  that  no  such 
shocks  occur,  except  in  the  limit,  v  ->  »  .  However,  for 

p 

large  values  of  the  parameter  ~~  defining  the  surface  load, 
the  solutions  come  extremely  close  to  the  singular  values 
representing  a  shock,  so  that  computing  difficulties  occur. 

d)  Asymptotic  Solutions  near  Singularities. 

As  stated  in  the  previous  paragraph,  numerical  difficulties 
in  the  vicinity  of  0  =  <£  will  make  the  procedure  for  inte¬ 
gration  of  Eqs .  (2-45)  outlined  in  subsection  b  unsuitable  and 

inapplicable  if  the  values  of  3  and  y  are  sufficiently  close 
to  those  for  a  plastic  front,  3  =  3,  7  =  .  To  establish 

the  behavior  of  the  solution  of  Eqs.  (2-45)  in  such  cases,  let 

T  =  f  +  r| 

3  =  3  +  A  (3-34) 

0  =  0  +  £ 

where  r),  A  and  €  are  small  quantities.  Introducing  these 
expressions  into  the  determinantal  equation  and  into  the 
differential  equation  (20)  for  0'  ,  retaining  only  the  leading 

terms  in  each  of  the  new  variables,  one  obtains  for  v  ^  ■g 
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A2  -  a  T)2  =  -  a  e 

1  2 

(3-35) 

f»  =  i  +  L 

(3-36) 

Combination  of  Eqs .  (17),  (18)  and  (19)  permits  the  formulation 

of  a  relation  for  ~  =  a  1 

de  H  J 

de  =  -  A>  L  (3-37> 


These  three  equations 
three  unknowns  tj  ,  A 


dJ^ 

de 


govern  the  solution  in  terms  of  the 
and  L.  The  derivative  of  J  becomes 


_ 6(l+v)  kL _ 

(i-2v)-\|3+(3+I)* 


(3-38) 


while  s  and  s  can  be  found  from  the  yield  condition  and 
1  2 

from  Eq.  (34)  once  A  is  known.  A  solution  of  the  above  equa¬ 
tions  is  valid  only  if  the  inequality  (28) ,  L  >  0  ,  is  satis¬ 
fied. 


The  coefficients  are 

_  288  (l+v). 
1  "  l-bv 


H 


_  48(l+v) 
l-2v 


3(l-v)V2  _  -j 


(l+v) 


(3-39) 
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} 


=  115. 

l- 

_  g( l+v 
l-'Hv 


ST 


The  terms  (l-8v)  in  the  denominator  of  several  coefficients 

necessitate  exclusion  of  the  case  v  =  1/8  . 

The  three  asymptotic  equations  (35-37)  are  nonlinear  and 

*) 

do  not  have  closed  solutions  .  Numerical  integration  near 
e  =  A  =  r]  =  0  would  again  encounter  difficulties  because  of 
the  presence  of  the  same  type  of  singularity  which  occurs  in 
the  original  differential  equations.  However,  Appendix  A  shows 
that  simple  expressions  exist  which  describe  the  asymptotic 
behavior  of  the  solution  of  Eqs.  (35-37)  near  e  =  A  =  t]  =  0. 
These  equations  being  asymptotic  approximations  for  the 
original  differential  equations,  the  behavior  of  the  solutions 


*) 

By  elimination  of  L  and  A  one  can  find,  as  alternative 
to  Eqs.  (35)- (37),  a  first  order  differential  equation 


c  n  -  e 

— -  (3-40) 

c  r\2  -  e 
2 


where 


c 

2 


ai(1+a4) 

a 

2 


(3-41) 


valid  provided 


a  /  -  1  (i.e.  v  /  1/2). 
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of  the  latter  is  also  described  by  the  results  obtained  in 
Appendix  A.  It  is  shown  that  Eqs .  (35-37)  for  v  /  1/8  have 
limiting  solutions  consisting  of  separate  branches.  Because 
the  signs  of  several  of  the  coefficients  a^  differ  for  v  ^  1/8, 
the  two  cases  have  different  character. 

1/2  )  v_>  1/8-  In  this  case  solutions  exist  only  for  e  /  0, 

which  is  due  to  the  charact'  r  of  Eq.  (35)  .  Noting 

a  /  0  ,  a  )>  0,  the  left  hand  side  of  the  latter  is  a  sum  of 
1  2 

squares  and  therefore  positive  and  no  real  solutions  A,  rj  can 
exist  for  e  )>  0. 

Two  types  of  limiting  solutions  exist,  given  by  Eqs. 
(A-3,4,7)  and  Eqs.  (A-8,9,14),  respectively. 

In  solutions  of  Type  1,  rj  is  proportional  to  7  -  e  , 
while  A  ,  which  is  proportional  to  e  ,  is  much  smaller  than  t), 
so  that  approximately  A  -  0  .  The  reverse  applies  for  solu¬ 
tions  of  Type  2,  where  A  is  proportional  to  J  -  e  and 
h  -  0  .  The  sign  of  the  square  roots  is  arbitrary,  so  that  in 
each  case  the  solution  for  the  respective  quantity  has  two 
branches . 

Whenever  the  numerical  integration  of  the  original  differ¬ 
ential  equations,  (2-45),  approaches  0  =  0,  the  values  r) 
and  A  must  approach  one  of  the  asymptotic  solutions.  In 
Section  5a,  where  numerical  cases  are  discussed,  the  actual 
and  the  asymptotic  solutions  for  A  are  shown  in  Fig.  15  for 
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1 


a  typical  case  where  v  /  1/8  .  The  values  of  the  stress 

deviators  being  defined  by  3  =  3  +  A  ,  and  A  being  a 

small  quantity,  s  and  s  at  the  terminal  point  of  the 

1  2 

integration  (i.e.  at  the  end  of  the  plastic  region)  must  have 
a  value  close  to  that  for  3=3,  viz. 


s 

l 


+  2-ALJs 


(3-42) 


The  signs  in  these  expressions  depend  on  the  branch  used  in 
aPProaching  0^0  .  The  value  of  j  is  to  be  determined  from 

l 

Eq.  (38).  Regardless  of  the  type  of  solution,  Eqs .  (A-4)  and 

(A-9)  indicate  that  L,  and  therefore  J'  ,  are  proportional 

l 

to  —  .  As  the  terminal  point  of  the  integration  moves  closer 

to  e  =  0,  the  integral  of  J'  ,  i.e.  the  value  j  ,  increases 

l  l 

without  upper  bound.  While  ,  0  and  7  are  practically 

constant  near  the  singularity/  a  very  large  change  in  j 

1 

occurs  in  a  very  small  angular  region.  The  magnitude  of  this 

change  depends  on  the  stopping  point,  i.e.  on  the  end  of  the 

plastic  region.  In  a  plastic  shock  J  is  discontinuous,  while 

1 

s.  .0  and  7  remain  constant  and  the  solution  described  above 
1 

is  quite  similar  to  such  a  shock,  except  that  the  change  in 
J  occurs  in  a  small,  but  finite  region  of  0  .  The  solution 

-L 

describes,  therefore,  a  "plastic  shock"  of  finite  thickness. 
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0  <(  v  1/8  .  In  this  case  a  /  0  ,  so  that  Eq.  (35)  permits 

solutions  without  restriction  on  e  .  A  plastic  region  may 
therefore  contain  the  point  0  =  0.  It  is  shown  in  Appendix 
A  that  for  e  /  0  solutions  of  Type  1  apply,  Eqs .  (A-3,4,7), 

where  T]  is  proportional  to  sfe  ,  while  A  <(  <(  r)  is  linear  in 
e  .  For  e  0  the  solution  is  of  Type  2,  A  being  propor¬ 
tional  to  \l  -  e  ,  while  T)  /  <(  A  is  linear  in  £  .  A  solution 
which  extends  continuously  to  both  sides  of  0  =  0  must  there¬ 
fore  follow  first  the  one,  then  the  other  type  of  limiting 
solution.  However,  the  two  are  not  continuous  in  the  deriva¬ 
tives,  so  that  an  actual  solution  cannot  follow  either  one 
asymptotically  to  the  origin  e  =  0  as  was  the  case  for  v  /  1/8, 
An  actual  solution  must  have  a  smooth  transition,  contained 
in  Eqs.  (35-37),  but  lost  in  Appendix  A  due  to  the  approxi¬ 
mations  required.  Figures  17  and  18,  presented  in  Section  5a 
in  connection  with  typical  numerical  results,  show  the  actual 
and  the  asymptotic  solutions. 

The  details  of  the  transition  could  be  obtained  by  numer¬ 
ical  integration  of  Eqs.  (35-37).  Even  without  such  inte¬ 
gration,  one  can  already  state  qualitatively  that,  just  as  for 
v  y  1/8,  the  quantities  s^,  0,  y  will  undergo  only  minor 


changes  near  0=0 


:he  — ■*” 


v  aiuc  vj 


V\  /"<  *  .  »  t  v*  ».»i  1  "1  y— < 

,  i iv-/vVc:  v  ci  /  vVxxx  i xcii l'-j G 


radically,  it  becomes  larger  without  bound,  when  the  solution 
passes  closer  to  the  point  rj  =  A  =  e  =  0  .  Again,  the  actual 
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solution  in  the  vicinity  of  the  singularity  may  be  said  to  be 
a  plastic  shock  of  finite  thickness. 

When  the  asymptotic  solutions  were  obtained  it  was  ex¬ 
pected  that  the  numerical  integration  of  the  original  differ¬ 
ential  equations  could  be  carried  sufficiently  close  to  the 
singularity,  so  that  the  range  of  validity  of  the  two  solutions 
overlaps.  Due  to  the  very  severe  singularity  this  expectation 
was  not  born  out  by  the  facts.  To  obtain  a  satisfactory  overlap 
of  solutions,  approximate  equations,  similar  to  Eqs.  (35-37), 
but  retaining  higher  order  terms,  were  formed,  and  integrated 
numerically.  The  details  are  given  in  [8]  •  While  the  solutions 
in  closed  form  obtained  in  Appendix  A  give  a  good  qualitative 
understanding  of  the  shock  front  of  finite  thickness  described 
above,  they  are  not  sufficient  to  find  quantitative  results, 


which  can  be  obtained  from  the  analysis  in  [8] 


For  this 


reason,  combined  with  the  complexity  of  the  refined  equations 
required  for  the  special  case  v  =  1/8  ,  no  attempt  was  made  to 
find  an  asymptotic  solution  for  this  case  in  closed  form. 

An  important  conclusion  which  may  be  drawn  from  the  char¬ 
acter  of  the  asymptotic  solutions  concerns  the  question  whether 
a  plastic  shock  may  occur  in  the  interior,  or  at  the  boundary 
of  a  plastic  region.  If  this  would  occur  the  change  in  value 

of  J  in  an  interval  including  the  shock  location  0  would  be 

1  dJi  .  . 

infinite,  because  the  integral  of  is  divergent  at  the 
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** p-=- 


point  A  =  T)  =  e  =  0  .  This  of  course  can  not  occur  in  an 

P0 

actual  case  where  the  surface  load  -j—  is  finite,  so  that  only 
the  possibility  of  a  shock  front  between  nondissipat ive  regions 
is  to  be  considered  in  the  construction  of  solutions. 


4. 


CONSTRUCTION  OF  SOLUTIONS. 


In  Section  3  a  number  of  partial  solutions  were  obtained 
from  which  the  solution  of  the  complete  boundary  value  problem 
is  now  to  be  constructed.  Section  3b  gives  the  differential 
equations  for  the  determination  of  the  stresses  and  velocities 
in  plastic  regions;  Section  3a  indicates  that  all  unknowns 
in  nondissipat ive  regions  are  constants,  except  for  discon¬ 
tinuities  of  a  prescribed  nature  at  the  locations  and  0p  . 

In  addition,  as  discussed  in  Section  3c,  a  shock  front  with 
plastic  deformation  may  occur  at  the  location  0. 

In  Section  2  boundary  conditions,  and  additional  require¬ 
ments,  which  the  solution  must  satisfy,  were  formulated  and 
discussed.  Equations  (2-16  and  17)  for  the  prescribed  surface 
load  in  terms  of  the  variables  s^  ,  J  ,  y  and  0  require 
either 


J.  (tt) 

s ^ ( tt)  +  — 3 —  =  -  PQ  ,  0(tt)  =  2  (4-la) 


,  x  Ji(*) 

s2 +  — ^ =  “  Po  *  6(tt)  =  0  (4-lb) 

A  further  boundary  condition  requires  that  all  quantities 
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must  vanish  for  0^0 p  ,  Eq.  (2-18) .  This  condition,  in 
conjunction  with  the  fact  that  a  plastic  region  or  a  plastic 
shock  can  exist  only  in  locations  where  the  yield  condition 
is  satisfied,  permits  the  conclusion  that  the  change  in  stress 
from  vanishing  values  for  0  <(  0p  to  nonvanishing  values 
must  be  ncndissipative.  However,  in  nondissipative  regions 
the  stresses  are  constant,  except  for  discontinuities  at 
0=0 p  or  0=0  .  A  solution  in  which  plastic  deformations 
occur  at  all  can  therefore  siart  only  in  one  of  the  two  ways 
described  below. 


Case  1.  Discontinuities  occur  at  the  P-  and  S-front 
where  the  discontinuity  at  0p  satisfies  the 
inequality 


a 

l 


< 


y/~3(  1- v)  k 
l-2v 


(4-2) 


while  the  discontinuity  At  at  0g  is  of  such 
magnitude  that  the  yield  condition 


F 


0 


(4-3) 


is  satisfied. 


The  symbols 


(-) 


indicate 


approach  from  above  or  below,  respectively. 


Case  2 .  A  discontinuity  occurs  at 
cribed  by  Eqs .  (3-7,8),  so  that  the 


is  satisfied  for  0  =  0 


=  *(+) 


i  .  e . 


0 


(  +  ) 


vA(  l-v)k 
l-2v 


the  P-front,  des- 
yield  condition 


(4-4) 


In  Case  1,  plasticity  may  occur  only  in  locations  0  )>  0  , 

b 

while,  in  Case  2  it  may  occur  already  for  0  <(  0  <(  0 

P  S 

As  a  next  step  in  the  search  for  solutions  it  is  helpful 
to  consider  the  latter  as  function  of  the  nondimensional  sur- 

Po 

face  pressure  ,  while  Poisson's  ratio  v  and  the  velocity  V 

P0 

are  considered  constant.  For  sufficiently  small  values  of 

po 

the  solution  must  be  entirely  elastic,  but  as  increases  plastic 

regions  must  occur  and  should  form  a  gradually  changing  pattern. 

PE 

Based  on  [1]  one  can  find  the  limiting  value  ,  above  which 
entirely  elastic  ■•■•.olutions  are  no  longer  possible. 

The  elastic  solution,  shown  in  Fig.  7,  has  two  discon¬ 
tinuities  at  0 p  and  0g  with  regions  of  constant  stress 
between  0p  and  0g  and  between  and  the  loaded  surface 

0  =  7f  .  Depending  on  the  values  of  the  parameters  v  and  ~ 

CP 

the  yield  condition  may  be  reached  in  either  of  the  two  regions 

pp 

resulting  in  different  expressions  for  ~  .  If  the  region 
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0  )>  0O  controls , 


_E 

k 


M 


2 


3N2  -  3N  cos  20  _  +  (l-v+v2)  cos2  20_ 

o  b 


(4-5a) 


where 


N  =  •£ 


cos  20c  +  (l-2v)  cos  2(0q  - 


♦r) 


(4-5b) 


while 


_ _ ]2 

( l-2v) 2  cos2  20 s 


(4-6) 


if  the  region  0  \  0g  controls.  The  decision  which  region 
controls  can  be  made  by  comparing  the  values  given  by  F.qs .  (5) 
^.id  (6),  the  smalle..  one  controlling.  Designating  as  Range  I 
the  combination  of  values  v  and  —  where  Eq.  (5a) controls , 

Cp 

one  finds  that  in  this  range 


( l-2v)  2 
(l-v)"(  l-3v) 


(4-7) 


-39- 


') 

i 

•4 


The  remainder  of  the  range  will  be  designated  as  Range  II. 

PE 

Both  ranges  are  shown  in  Fig.  8.  The  limiting  values  —  are 


shown  in  Fig.  9  for  several  values  of  v  as  function  of  . 

pE  cp 

If  the  surface  load  exceeds  the  value  —  by  a  sufficiently 
small  amount  the  elastic-plastic  solution  should  differ  only 
slightly  from  the  elastic  one,  which  can  be  used  to  predict 
the  character  of  the  solution.  Because  the  situations  differ, 
the  Ranges  I  and  II  must  be  discussed  separately. 


a)  Solutions  in  Range  I. 

PE 

In  this  case  is  given  by  Eq.  (5a)  and  the  yield  stress 

in  the  elastic  solution  is  reached  only  in  the  region  0  )>  0g  . 
The  discontinuity  a  satisfies  then  the  inequality  (2)  and 
continuity  requires  that  this  inequality  will  still  apply  for 


E 


a  range  of  values  )>  ,  where  pL  is  a  limiting 

value,  not  yet  known.  In  this  range  the  start  of  the  solutions 
will  be  according  to  Case  1  and  plasticity  can  therefore  occur 
only  in  the  region  0  )  0g 

Using  an  indirect  approach,  the  determination  of  plastic 

PQ  /  PL 

solutions  for  values  \  begins  with  the  selection  of  a 
pair  of  starting  values  a  and  At  near  those  for  the  limit¬ 
ing  elastic  case.  Experience  and  continuity  considerations 
indicate  f.hat  the  value  \  a  |  should  be  larger  than  |  a  | 
corresponding  to  p  given  by  Eq.  (5a).  A  plastic  region  can 
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start  only  at  a  point  0  which  is  a  root  of  Eq.  (3-13). 

Inequalities  derived  in  Appendix  B,  Eq.  (B-23)  indicate  that 

this  equation  has  one,  but  only  one  root  0  )>  0_  .  Starting 

1  S 

integration  at  0  =  0  the  solution  in  the  interior  of  the 

plastic  region  is  determined  from  Eqs .  (3-18  to  27).  The 

plastic  region  can  be  extended  as  long  as  Eq.  (3-27)  gives 

values  L  )>  0,  but  the  plastic  region  may  be  terminated  at 

will  at  any  earlier  location  0  .  The  solution  for  0  )>  0 

2  2 

is  then  nondissipative ,  i.e.  all  quantities  are  constant.  If, 

therefore,  in  the  process  of  forward  integration,  a  value  9  is 

encountered  which  satisfies  Eq.  (1) ,  the  plastic  region  is 

terminated  and  a  solution  for  one  value  of  the  surface  pressure 
P0 

has  been  obtained.  Repeating  this  process  with  gradually 

increasing  starting  values  |  o  |  the  whole  spectrum  of  values 

satisfying  the  inequality  (2)  can  be  explored.  Solutions,  if 

any,  obtained  in  this  manner  will  have  the  configuration  shown 

in  Fig.  10,  i.e.  discontinuities  at  0  and  0  and  a  plastic 

P  ^ 

region  7T  )>  0  >  0  >  0  •  There  will  be  an  elastic  region 

2  l  S 

of  constant  stress  from  0  to  0  and  two  neutral  regions  to 

P  O 

either  side  of  the  plastic  one. 

Postponing  the  discussion  of  uniqueness  and  of  alternative 

configurations,  one  can  proceed  in  a  similar  manner  when  the 

solution  begins  at  0  =  0p  as  indicated  in  Case  2.  Starting 

with  a  value  a  according  to  Eq.  (3) ,  the  yield  condition  is 

l 
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,  so  that  the  determinantal 


satisfied  for  any  value  0  )>  0^ 

equation  (3-13)  according  to  Eq.  (B-23)  now  has  two  roots, 

0  ^  0g  ,  at  which  plastic  regions  may  start.  Both  roots 

must  be  explored.  If  the  larger  one,  0  )>  0  ,  _eads  to  a 

l  S 

solution,  it  has  a  configuration  as  shown  in  Fig.  11a.  Start¬ 
ing,  alternatively,  with  the  smalle1  oot ,  0D  <  0  <  0C  ' 

Pi  S 

several  possibilities  are  to  be  investigated.  The  integration 

may  be  continued  as  long  as  l  )  0  to  see  if  a  value  0=0  or 

—  can  be  reached.  The  configuration  of  such  a  solution,  if 

any,  is  shown  in  Fig.  11b.  Alternatively,  the  plastic  region 

can  be  terminated  at  will  at  a  point  0  <(  0  where  0  ^  0,  -S 

2  S  d 

The  plastic  region  will  then  be  followed  by  a  neutral  one  for 
values  0  )>  0  .  The  inequalities  on  the  roots  of  Eq.  (3-13) 


indicate  that  there  is  just  one  more  root  0  )>  0  ,  when  a 

3  S 

second  plastic  region  can  begin.  Starting  integration  at 
this  point  may  lead  to  a  terminal  location  0  ,  where  0=0 

4 

or  .  The  configuration  of  such  a  solution,  if  any,  Fig.  11c, 

contains  a  p-front  and  two  plastic  regions,  separated  by  three 

neutral  regions  There  are,  however,  further  possibilities. 

The  neutral  region  0  )>  0  which  follows  the  first  plastic 

2 

one  may  be  terminated  at  0S  by  an  elastic  change  in  shear, 

At  -  which  is  restricted  in  sign  and  intensity  by  the  yield 


condition.  If  At  is  such  that  F 


0, 


(  +  ) 


0  ,  the  region 


rr 


0  y  0^  becomes  elastic.  This  might  permit  values  6  =  0,  75 
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at  the  surface,  the  corresponding  solution  having  the  config¬ 


uration  of  Fig.  12.  Finally,  the  important  case  must  be  con¬ 


sidered  where  the  value  of  At  is  such  that  F 


0 


(+) 


vanishes 


again,  a  situtation  discussed  in  Section  3-a-2.  In  the  latter 

case  there  is  again  a  neutral  region  for  0  )>  0  ,  which  can 

b 

be  followed  by  a  plastic  region  because  Eq.  (3-13)  has  a  root 
0  >  0O  giving  a  starting  location.  The  configuration  of  a 

3  0 

solution  obtained  in  this  manner  is  shown  in  Fig.  13. 


b)  Solutions  in  Range  II, 


•  T? 

In  this  range  is  given  by  Eq.  (6)  so  that  in  the 

,  .  .  .  Pn  Pp 

limiting  case  — ~  yield  is  just  reached  in  the  region 

K.  iC 

0p  <  0  <  0S 


The  discontinuity  a  at  the  P-front  must 

therefore  satisfy  Eq.  (4) ,  which  will  also  hold  for  neigh- 

P~  PE 

k 


boring  elastic-plastic  solutions  where  —  exceeds  slightly. 


These  solutions  will  therefore  start  at  0  =  0p  according  to 

PQ  PEP 

Case  2.  In  the  limiting  solution  for  the  region 

0  >  0S  is  below  yield  and  continuity  requires  this  to  hold 

in  neighboring  elastic-plastic  solutions,  so  that  the  plastic 

region  must  lie  .in  the  range  0  <(  0  0  <(  0  .  The  con- 

p  i  2  S  „ 

struction  of  solutions,  Fig.  12,  begins  exactly  as  for  -  -  )>  — J- 

K  Jc 

For  each  terminal  point  0£  the  strength  of  the  discontinuity 
At  at  the  shear  front  is  determined  by  the  requirement  that 


0=0  or  ^  subject  to  the  limitation  F 


0 


(+) 


<  0 


When 
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the  required  value  At  violates  this  condition  a  second  plastic 
region  for  0  )>  0g  is  needed,  i.e.  the  configurations  shown 
in  Figs.  11c  and  13  are  to  be  investigated. 

c)  Alternative  Solutions,  and  Considerations  of  Uniqueness 
and  Existence. 

In  the  absence  of  a  uniqueness  theorem  it  is  vital  to 
demonstrate  that  configurations  other  than  those  shown  in 
Figs.  10  to  13  can  not  lead  to  solutions.  It  has  been  shown 
in  Appendix  B  that  Eq.  (3-13) ,  which  is  also  the  equation  for 
the  characteristic  directions,  has  for  a  given  state  of  stress 
one  root  0  ,  no  more  no  less,  in  each  of  the  intervals 
0p  <\  0  <\  0g  and  0g  <(  0  <(  n  .  If  a  plastic  region  ends  at  a 
location  0..  in  one  of  these  intervals,  the  state  of  stress 
in  the  remainder  of  the  interval  0^0^  is  necessarily  neutral 
and  uniform,  and  equal  to  the  one  at  the  terminal  point  0^  of 
the  plastic  region.  0^  is  therefore  the  only  solution  of 
Eq.  (3-13)  for  this  state  of  stress  in  the  particular  interval, 
and  no  more  than  one  plastic  zone  can  therefore  occur  in  any 
interval . 

In  Section  3-c  the  possibility  of  discontinuous  plastic 
shock  fronts  has  been  indicated  and  their  occurrence  must  be 
considered.  First,  it  has  been  concluded  at  the  end  of  Section 
3-d  that  a  plastic  shock  can  not  occur  in  the  interior  of  a 
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plastic  region  nor  at  its  end,  so  that  a  plastic  shock,  if  it 
occurs,  must  lie  between  nondissipative  regions  and  be  quite 
separate  from  a  continuous  plastic  region.  Further,  such  fronts 
can  not  occur  in  an  interval  0  <(  0  \  0  or  0  /  0  <(  rr  where 

Jr  o  o 

a  plastic  zone  occurs,  because  Eq.  (3-13)  which  is  satisfied 
in  the  location  of  a  plastic  shock  would  then  have  two  roots 
in  the  same  interval.  This  reasoning  leaves  only  the  possi¬ 
bility  of  configurations  similar  to  Figs.  10  to  13,  where  one 

of  the  plastic  zones  is  replaced  by  a  plastic  shock.  For  v  )>  1/8, 

-  .  *) 
where  0  \  0  ,  this  can  not  occur  (except  in  the  limit  v  ->  <») 

b 

because  Eq.  (3-31)  requires  that  the  principal  stress  be  normal 
to  the  shock  front  0  ,  while  the  actual  principal  stress  in  the 

neutral  region  behind  the  P-front  is  normal  to  the  latter. 

For  v  <(  1/8  ,  0  )>  0  ,  so  that  configurations  similar  to  Figs. 

O 

10  to  13  might  occur  where  the  plastic  region  degenerates  into 

a  shock  at  0  )>  0  .  This  can  not  happen,  however,  because 

b 

the  direction  of  the  principal  stress  behind  the  shock  front, 
at  0  )  0^+^  ,  would  have  to  be  normal  to  this  front,  which 

contradicts  the  requirement  0=0  or  ~  ,  on  the  surface, 
except  in  the  limiting  case  v  =  c »  •  Discontinuous  shock 
fronts  can  therefore  not  occur  at  all  for  finite  values  of  the 
velocity  V,  but  values  of  y  and  (3  where  the  conditions  (3-31) 

*) 

In  this  case  the  surface  pressure  is  applied  simultane¬ 
ously  everywhere  on  the  surface,  producing  the  trivial  result 
of  a  P-wave  followed  by  a  plastic  shock,  both  having  horizon¬ 
tal  plane  fronts. 
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are  nearly  satisfied  are  encountered.  The  asymptotic  behavior 
of  the  solution  near  0  =  0  in  such  cases  was  studied  in 
Section  3-d,  and  examples  are  given  in  Section  5. 

The  preceding  discussion  shows  that  for  finite  values  of 
V  no  plastic  shock  front  can  occur  and  that  there  can  be  no 
more  than  one  plastic  region  in  each  of  the  intervals 
0p  \  0  <(  0S  ,  0g  {  0  (  7r  .  Combined  with  elastic  disconti¬ 
nuities  at  the  P-  and  S-fronts,  only  the  limited  number  of 
configurations  shown  in  Figs.  10  to  13  are  possible. 

The  numerical  analysis  by  digital  computer  was  set  up  to 
investigate  all  possible  alternatives,  i.e.  the  configuration 
according  to  Fig.  10  if  the  starting  value  a  satisfies  Eq. 

l 

(2)  and  any  of  the  alternatives  shown  in  Figs.  11  to  13  if  a 
satisfies  Eq.  (3) .  While  none  of  the  configurations  shown  in 
Figs.  11  a-c  ever  furnished  a  solution,  no  general  proof  per¬ 
mitting  elimination  of  these  cases  is  available. 

In  Range  I  solutions  which  start  according  to  Case  1, 
have  the  configuration  of  Fig.  10.  For  fixed  values  of  v  and 
V,  these  solutions  form  a  family  which  depends  on  one  parameter, 

the  selected  starting  value  |  c  |  )>  |  |  .  It  was  found  that 

Po  , 

the  surface  load  increases  monotonically  with  |  o  |  until 

the  limit,  Eq.  (4)  for  o  is  reached,  which  leads  to  a  limit- 

1  PL 

ing  value  of  the  surface  load .  However,  no  analytical 

Po 

proof  of  the  monotonic  increase  of  is  available. 
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The  solutions  found  for  Range  I,  which  start  according 


to  Case  2,  had  always  the  configuration  shown  in  Fig.  13. 

These  solutions  also  depend  on  one  parameter,  viz.  the  stopping 

point  0  of  the  plastic  region  between  0_  and  0„  If  0 

2  PS  2 

is  selected  only  slightly  larger  than  0  ,  the  solution  must 

obviously  be  very  close  to  the  limiting  one  for  Case  1,  so 

P0  ~  PL 

that  in  such  a  case  >  —  and  there  is  a  smooth  transition 

from  the  configuration  according  to  Fig.  10  to  that  of  Fig.  13. 

Po 

The  numerical  analysis  indicated  that  the  surface  load 

increases  monotonically  with  0  .  As  0  approaches  a  limit- 

2  2  P 

ing  value  the  surface  load  goes  to  the  limit  ->  »  ,  for 
reasons  explained  in  Section  3-d. 

In  Range  II,  only  solutions  which  start  according  to  Case 
2  were  found,  their  configurations  being  as  shown  in  Figs.  12 

pE  .  p0  pL  pr 

and  13.  Figure  12  applied  for  values  \  where 

is  a  bound.  The  corresponding  family  of  solutions  depends  on 

Pt 

the  stopping  point  0  of  the  plastic  region.  The  bound  -~ 

2  K 

is  reached  when  the  elastic  region  for  0  )>  0  becomes  neutral. 

po 

For  larger  values  of  Fig.  13  applies  and  all  statements 

made  in  Range  I  for  this  case  apply. 

In  Range  I  as  well  as  in  Range  II,  combination  of  all 
solutions  obtained  numerically  furnished  one,  and  only  one 

pO  pe 

solution  for  each  value  of  )>  .  However,  no  general 

proof  is  available  that  this  must  oe  so.  Existence  and 
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uniqueness  of  the  solutions  obtained  must  therefore  be  demon¬ 
strated  for  each  combination  of  values  v  and  —  by  actual 

CP 

computation  of  the  families  of  solution  according  to  the  con¬ 
figurations  shown  in  Figs.  10  to  13. 
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5. 


RESULTS  AND  CONCLUSIONS. 


a)  Discussion  of  Typical  Numerical  Results 

The  numerical  integration  of  the  simultaneous  differential 

equations  (3-18  to  21)  in  plastic  regions  was  accomplished  by 

a  Runge-Kutta  forward  integration  scheme  of  fourth  order  [9] 

Computations  were  programmed  in  Fortran  II  for  an  IBM  7090 

digital  computer.  Only  typical  results  representing  the  various 

configurations  will  be  given  here  for  the  case  of  the  velocity 

PQ 

V  =  1.25  Cp  for  several  values  of  —  and  v  •  An  extensive 
tabulation  of  additional  results  is  given  elsewhere,  [8]. 
Significant  differences  :n  the  solution  for  large  values 
occur,  depending  on  whether  v  ^  1/8  ,  so  that  examples  for 
both  situations  are  presented. 

p 

For  v  =  1/4,  V  =  1.25  c  the  limiting  value  —■  up  to 

lr  ,K 

•  Pp 

which  the  solution  is  entirely  elastic  is  —  =  2.50.  Figure 

K 

14a  shows  the  detail  of  a  solution  for  a  slightly  higher  value 

P0  .  Pp 

of  the  pressure,  —  =  2.61  /  — ^  .  The  configuration  has  one 

plastic  region  in  agreement  with  Fig.  10.  The  limiting  value 

Pt 

for  this  configuration  is  ~  =  2.65  •  When  the  surface  pressure 

exceeds  this  value,  the  configuration  contains  two  plastic 

zones  and  is  of  the  type  shown  in  Fig.  13.  The  details  for  the 
o 

*  /*» 

value—  =  3.58  are  shown  in  Fig.  14b.  Further  solutions  for 

JC 

p0 

other  values  of  can  be  obtained  by  varying  the  end  point 
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0  of  the  lower  plastic  region. 

However,  computational  difficulties  arise  when  the  end 
point  0^  approaches  the  value  0  ,  which  has  in  the  present 
example  the  value  0  =  143.3957°  .  As  explained  in  Section 
3-d  a  very  rapid  change  occurs  in  the  quantity  j  ,  nearly 

l 

a  shock  front,  when  0  ->  0  ,  with  the  result  that  the  corres- 

2 

ponding  surface  pressure  will  be  large.  Figure  15  illustrates 
the  result  of  the  measures  taken  to  obtain  numerical  results 
in  this  range.  First  the  program  for  the  numerical  integration 
was  revised  for  double  accuracy  (16  digits) .  When  it  was  found 
that  this  was  not  adequate  to  continue  solutions  sufficiently 
close  to  0=0  to  obtain  agreement  with  the  asymptotic  solu¬ 
tions,  expansions  including  higher  order  terms  than  used  in 

*) 

Appendix  A  were  made  and  the  resulting  approximate  differen- 

** ) 

tial  equations  were  integrated  numerically  .  This  procedure 
was  successful  as  illustrated  by  Fig.  15.  It  was  stated  in 
Section  3-d  that  for  v  >  1/8  there  are  two  types  of  asymptotic 
solutions  and  the  actual  solution  may  approach  either  one. 

Figure  15  shows  the  quantity  A  =  (3  -  3  as  function  of  e  =  0  -  0. 
The  solution  of  the  approximate  differential  equations  approaches 

*) 

The  details  of  these  numerical  computations  are  contained 
in  [8]  . 

**)  _ 

The  expansions  indicate  that  near  0=0  cancellations  of 
leading  terms  up  to  high  order  occur,  which  explains  the  diffi¬ 
culty  encountered. 
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the  negative  branch  of  the  asymptotic  solution  of  Type  2  for 
very  small  values  of  e  .  The  curve  shown  is  part  of  the  solu- 

PQ 

tion  for  =  10, 5  •  The  integral  of  the  original  differential 

equations  even  with  double  precision  approaches  the  asymptotic 

solution  not  well  and  is  in  this  range  inaccurate.  However,  the 

solution  of  the  original  and  of  the  approximate  differential 

equations  agreed  very  well  for  e>  5  x  10~5,  which  fact  could  not 

be  shown  in  the  scale  of  Fig.  15 .  The  complete  results  for  the 
p 

case  -rr-  =  10.5,  corresponding  to  Fig.  15,  are  shown  in  Fig.  l4c . 
Solutions  for  other  values  of  v  and  V  are  essentially 

p 

similar,  except  for  solutions  for  large  values  of  ~  when 

K 

v  <(  1/8  ,  because  the  asymptotic  behavior  near  0  =  0  changes. 

To  illustrate  such  a  case,  Figs.  16a-c  give  solutions  for  v  =  0, 

p  p 

V  =  1.25  cp  ,  where  —  =  l.ij.1  and  =  2.40  .  Figures 

Pp  P  Pt 

16a,  b  apply  respectively  for  —  (  —  =  2. 17  <(  •—  and 

Pt  P  ^  p 

<(  ~  =  3.70  .  Figure  16c  finally  applies  for  —■  =  8.02 

a  case  where  the  numerical  integration  of  Eqs.  (3-35  to  38) 

encountered  computing  difficulties. 

As  discussed  in  general  in  Section  3-d,  in  the  last  case, 

the  singular  angle  0  is  situated  in  the  interior  of  the  upper 

plastic  region,  and  the  asymptotic  expressions  for  A  =  (3  -  3 

r|  =  7  —  tt/2  as  function  of  e  =  0  -  0  differ  for  e  <(  0  and  } 

e  y  0  ,  and  are  discontinuous  at  e  =  0  .  Figures  17,  18  show 

A  and  t)  given  by  the  asymptotic  expressions  and  the  actual 
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integral  obtained.  For  the  value  =  8.02  considered,  the 
approximate  differential  equations  obtained  by  expansion  of 
Eqs.  (3-35  to  38),  [8]  ,  had  to  be  used.  Figure  17  also  shows 

P0 

a  solution  for  a  lower  value  where  Eqs.  (3-18  to  21)  could 
still  be  integrated  using  double  accuracy.  As  stated  earlier, 

P0 

solutions  for  various  values  of  the  surface  pressure  are 

obtained  by  changing  the  end  point  0  of  the  lower  plastic 

2 

region.  In  the  range  under  discussion  the  computation  becomes 

very  sensitive  to  small  changes  in  0  .If  an  upper  bound, 

2 

which  can  only  be  found  numerically,  is  exceeded  the  integra¬ 
tion  no  longer  leads  to  a  solution  because  the  value  0  instead 
of  approaching  the  value  0  =  9O0  moves  away  from  it.  Figure 
17  shows  also  one  of  the  integrals  which  does  not  lead  to  a 
solution. 

Solutions  in  Range  II,  i.e.  for  values  of  v  and  — 

CP 

located  in  the  cross-hatched  area  of  Fig.  8,  do  not  differ  from 

p  pL 

those  in  Range  I,  except  when  <(  .  As  illustration  Fig. 

JC  K 

19  treats  the  case  v  =  .35,  V  =  1.25  cp  f°r  the  surface  pressure 

p  p  p 

<(  ~  =  4. 4l  <  =6.24  -  In  this  solution  an  elastic 

K  K  KL 

region  occurs  for  0  )>  0 

O 

b)  Conclusions . 

The  effects  of  a  step  pressure  pQ  progressing  with  con¬ 
stant  superseismic  velocity  V  )>  cp  on  the  surface  of  a  half- 
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space  have  been  determined  for  an  elastic-plastic  medium  obeying 
the  von  Kises  yield  condition.  The  steady-state  solutions 
obtained  satisfy  additional  conditions  selected  to  ensure  that 
the  solutions  are  asymptotic  ones  in  the  vicinity  of  the  front 
of  the  surf:.ce  load  for  transient  problems  of  the  type  shown 
in  Fig.  2. 

In  spite  of  the  lack  of  a  general  uniqueness  and  existence 

theorem,  a  unique  solution  was  obtained  for  each  combination 

.  V  Po 

of  the  significant  nondimens .ional  parameters  — —  ,  v  and  -rr 

P 

for  which  numerical  computations  were  actually  made.  Extensive 
numerical  results,  in  addition  to  the  typical  ones  discussed 
above  can  be  found  in  [8], 

The  elastic-plastic  configuration  of  the  solution  differs, 

there  being  three  cases,  Figs.  10,  12,  13,  depending  on  the 

P0 

value  of  the  nondimensional  parameters.  For  pressures 

PF  \r 

below  a  value  Tr=  ,  which  is  a  function  of  v  and  —  ,  the 

k  CP 

solutions  are  entirely  elastic.  For  larger  values  of  the 

PE  .  PQ  .  PL 

pressure,  in  a  range  \  —  ,  the  solutions  contain  one 

plastic  region,  the  configurations  being  shown  in  Figs.  10  and 

12.  The  former,  or  the  latter  applies  when  the  parameters  v 
V 

and  —  are  in  Range  I  or  Range  II,  respectively,  see  Fig.  8. 

CP  p 

For  values  of  the  nondimensional  pressure  above  a  limiting 

Pt  tt 

value  ~z —  ,  which  is  of  course  a  function  of  v  and  ”  , 

k  Cp 

two  plastic  regions  occur  as  shown  in  the  typical  configuration 
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Fig.  13. 


Two  significant  features  common  to  all  elastic-plastic 
solutions  found  should  be  noted. 

1.  While  in  one  dimensional  plane  or  spherical  problems, 
discontinuous  plastic  fronts  occur,  [ 5] f  [ 10] , 

they  do  not  exist  in  the  solutions  of  the  two  dimen¬ 
sional  problem  solved. 

2.  All  solutions  obtained  contain  an  elastic  discontinuity 
at  the  S-front,  in  addition  to  the  one  at  the  P-front. 
The  latter  occurs  also  in  one-dimensional  problems 

and  is  to  be  expected  from  general  considerations  of 
wave  propagation.  Such  considerations  also  permit 
prediction  of  the  possibility  of  a  discontinuity  in 
shear  similar  to  the  one  in  the  elastic  solution.  The 
solutions  found  show  that  this  discontinuity  occurs  in 
all  cases,  a  fact  which  can  not  be  established  by  purely 
qualitative  considerations. 

In  view  of  the  asymptotic  character  of  the  solutions  ob¬ 
tained,  it  must  be  expected  that  these  two  features  will  be 
found  also  in  steady-state  solutions  for  non-step  loads,  Fig. 

1,  and  in  transient  cases,  Fig.  2. 

The  present  paper  is  the  first  one  to  give  results  for 
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multidimensional  wave  propagation  in  elastic-plastic  media, 
except  for  purely  numerical,  finite  difference  schemes*)  in 
which  discontinuities  can  not  appear.  The  solutions  found  in 
the  present  paper  permit  checks  on  the  effectiveness  of  these 
numerical  schemes,  particularly  in  the  vicinity  of  disconti¬ 
nuities  in  the  actual  solution. 

As  a  by-product  of  the  steady-state  solution  for  the  step 
load,  the  character  of  the  partial  differential  equations  for 
the  general  case  was  examined  in  Appendix  B .  For  superseismic 
velocities  V > cp  the  equations  were  found  to  be  hyperbolic. 

The  method  used  in  this  paper  is  also  applicable  to  cases 
with  other  yield  conditions.  The  equivalent  problem  for  a 
yield  condition 


J 


2 


2  -r  2 


=  0 


(5-1) 


has  been  treated  concurrently  with  the  present  problem  [12], 


Computer  programs,  not  published,  or  given  limited  dis¬ 
tribution  only  have  been  developed  by  a  number  of  organizations 
for  various  purposes,  an  example  being  [11]  . 
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APPENDIX  A  -  Asymptotic  Solution  of  Equations  (3-35,36,37). 


The  asymptotic  equations  (3-35,36,37)  have  been  derived 
without  assumption  on  the  relative  magnitude  of  the  quantities 
i]  j  A  and  e  the  leading  terms  in  each  quantity  were  retained 
in  each  equation.  Further  simplifications  of  the  asymptotic 
equations  are  possible  by  studying  separately  the  various 
possibilities  for  the  relative  magnitude  of  -q  and  A 

If  rj  and  A  are  of  equal  magnitude,  Eq.  (3-35)  requires 
that  in  the  limit  e  •>  0  the  relations 

A  =  c  sTe  ri  =  c  s/~e~  (A-l) 

1  '  2 

hold,  where  c  and  c  are  constants,  This  assumption  leads, 
12 

however,  to  a  contradiction.  Substituting  Eq.  (A-l)  in  Eq. 
(3-37),  a  r|2  is  in  the  limit  negligible  compared  to  A  > 

3 

so  that 

L  =  -  h  (A-2) 

Substituting  this  expression  and  the  value  of  r]  in+o  Eq.  (3-36) , 

the  term  unity  is  negligible  compared  to  a  nL  •  The  result 

4  ' 

requires  a  =  -1  ,  valid  only  for  v  =  j/2  •  This  value 

4 

corresponds  to  an  incompressible  material  where  Cp  =  co  ,  so 
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that  the  superseismic  problem  to  be  studied  here,  does  not  exist. 
Solutions  according  to  Eqs.  (A-l)  do  not  apply  here. 

Solutions  of  Type  1.  If,  in  the  limit  e  ->  0,  T}  )>)>  A  ,  Eq. 
(3-35)  requires 


The  unity  term  being  small  compared  to 


1 

\fe 


,  Eq.  (3-36)  gives 


L 


1 


2a  e 

4 


(A-4) 


The  requirement  L  )> 0  restricts  the  sign  of  e  for  which  the 

solution  applies  to  sign  e  =  sign  a  .  Inspection  of  Eqs. 

4  a2e  \ 

(3-39)  shows  that  this  requirement  ensures  “  )  0  ,  so  that 

di 

Eq.  (3)  gives  real  values  of  T] 

Substitution  of  Eqs.  (A-3,4)  into  Eq.  (3-37)  gives  the 
differential  equation 


2a  e 

4 


dA 

de 


+  A 


a  a 
2,.  .3. 


(A- 5) 


the  general  solution  of  which  is 
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A  =  Ce 


(A-6) 


1 

2a 

4 


a  a  £ 

a  Tl  +2a  )' 

1  4 


where  C  is  an  arbitrary  constant.  However,  the  solution  (A-6) 
must  satisfy  the  premise  a/  <(  r)  where  t)  is  given  by  Eq.  (A-3) . 
If  the  constant  C  does  not  vanish  the  exponent  of  e  in  the 
above  expression  must  be  larger  than  ±/2  ,  or  —  /  1  .  Use 

H 

4 

of  the  expression  for  a  ,  Eq.  (3-39)  indicates  that  the  in- 

4 

equality  requires  v  )>  l/2  ,  an  impossible  requirement.  The 
constant  C  must  therefore  vanish,  and 


A 


a 

i 


a  ae 
2  3 

(1  +  2a  ) 
4 


(A- 7 ) 


a 

The  signs  of  a4  and  of  — 2-  which  govern  the  sign  of  e  ,  depend 

a  1 

on  the  value  of  v  •  The  solution  applies  for  e)>  0  when 
v  <(  1/8  ,  and  for  e  <(  0  when  v  /  1/8  .  The  denominator  in 
Eq.  (A-7)  does  not  vanish  in  the  range  0  V  <(  l/2  ,  so  that 

the  result  applies  except  for  the  previously  excluded  value 
v  =  1/8 


Solutions  of  Type  2.  If,  in  the  limit  e  ->  0,  A  //  T) 


Eq. 
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{ 


(3-35)  requires 


A  "  A/  V 


(A-8) 


a2  being  always  positive,  this  solutions  applies  only  for 
e  <(  0  .  Due  to  the  premise  A)  )  t]  the  term  a  t)2  in  Eq. 

3 

(3-37)  is,  in  the  limit,  small  compared  to  A  ,  giving 


(A-9) 


The  value  of  r|  is  to  be  determined  from  the  differential 
Equation  (3~36)  after  substitution  of  Eq.  (A-9) , 


£T|i  +  -q  =  e 


(A-10) 


If  a^  ^  -2  ,  the  general  solution  of  this  differential  equation 


Ce  2  + 


2  +  a 


(v  £  2/7)  (A-ll ) 


while  for  v  =  2/7  ,  when  a  =  -2 

4 


r)  =  Ce  +  eln  e 


(A-12) 


-59“ 


t 


where  C  indicates  an  arbitrary  constant. 

The  premise  A  /  )>  rj  limits  the  exponent  in  the  first  term 

a4 

of  Eq.  (A- 11)  ,  -  /  l/2  .  This  condition  is  not  satisfied 

for  v<l/8  ,  requiring  C  =  0  ,  but  it  is  satisfied  for  values 
of  v/l/8  .  Equation  (A-12)  for  v=  2/7  also  satisfies  the 
premise  A/  /  r\  ,  so  that  the  constant  C  does  not  vanish  in 
Eqs .  (A-11,12)  if  v  )>  l/8  .  For  v  <  i/8  the  solution  for  r) 

approaches  therefore  asymptotically  the  expression 


r]  -  2  +  a~  (v  <  i/8)  (A-13) 

4 

while  the  expressions  for  v  /  i/8  contain  an  arbitrary  constant 
and  are  not  fully  defined.  While  it  is  possible  to  find  a  range 
of  v  where  in  the  limit  e  ->  0  the  term  containing  C  is  negli¬ 
gible,  the  matter  need  not  ue  pursued  because  it  will  be  seen 
in  the  examples  that  for  v  /  i/8  it  is  sufficient  to  know  that 
r|  /  <(  A  ,  so  that  one  can  use  the  approximation 


A  >  >  r)  ~  0 


(v>  1/8)  (A- 14) 


in  conjunction  with  Eqs.  (A-8,9)  for  A  and  L 
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APPENDIX  B  -  Proof  for  the  Hyperbolic  Character  of  Equations 


(2-14,15)  and  Bounds  on  Characteristics. 

The  differential  equations  applicable  for  steady  state 
solutions  in  plastic  regions  are  the  six  relations  (2-14,15) 
supplemented  by  the  yield  condition,  F  =  0  .  To  prepare  the 
system  for  the  manipulations  required  to  obtain  the  charac¬ 
teristic  directions,  the  function  a  ,  the  derivative  of  which 

does  not  occur  in  Eqs .  (14),  (15),  may  be  considered  to  be  the 

c)A 

derivative  with  respect  to  £  of  another  function  A  ,  or  A.  = 

Further,  differentiation  of  the  yield  condition  converts  it 

^F 

into  a  differential  equation,  ^  =0  ,  a  total  of  seven  for 

seven  unknowns . 

These  seven  differential  equations  can  be  combined  with 
expressions  for  the  derivatives  of  the  unknowns  f  in  the 
direction  0  of  a  characteristic,  if  any, 

cos  0.  +  sin  0-  (B-l ) 

a£  i  dy  dr  v 

where  r  is  a  new  variable.  The  coefficients  of  the  derivatives 
with  respect  to  £  and  y  form  a  determinant  which  must  vanish 
in  a  hyperbolic  system  for  seven,  real,  characteristic  angles 
0^  .  It  is,  however,  not  necessary  to  derive  the  character¬ 

istic  equation,  because  the  seven  differential  equations  are 
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not  only  quasi-linear ,  but  also  homogeneous  in  the  derivatives 

d  f  S  f 

^  .  In  such  a  case  the  determinant  for  the  deter¬ 

mination  of  the  characteristic  angles  0^  is  necessarily  iden¬ 
tical  with  the  determinant  obtained  by  substitution  of  the 
relations  (2-24,25)  into  the  differential  equations.  This 
substitution  was  made  in  Section  2  after  a  change  of  variables 
which  cannot  be  material,  so  that  the  determinant  of  Eqs .  (2-45) 
must  vanish  for  characteristic  values  0^  .  This  is  precisely 

the  condition  which  leads  to  the  determinantal  Equation  (3-13) > 
which  applies  thus  for  the  values  0^  ,  too.  However,  it  will 

be  seen  that  Eq.  (3-13)  gives  four,  not  seven  roots  0^  ,  which 

difference  can  be  traced  to  disregarded  factors  containing 
powers  of  sin  0  ,  which  were  of  no  consequence  in  Section  3. 

That  there  is  a  triple  root  0^=0  ,  as  required  for  a  total 

of  seven,  may  be  proved  by  showing  that  the  matrix  of  the 

S  f 

coefficients  of  the  derivatives  ^  in  the  combined  Eqs. 

SF 

(2-14,15)  and  ^  =0  is  of  order  seven  but  only  of  rank  four. 

To  demonstrate  the  hyperbolic  character  of  the  differ¬ 
ential  equations  for  the  steady-state  the  roots  of  Eq.  (3-13) 
must  be  studied.  This  equation,  omitting  the  unessential 

factor  (s  +  s  )2  ,  may  be  written 

12 


C=[£cos  2y  +  (l-2v)(l-2X)]2+2[l+(i-2v)(l-2X)][(l+vXl-2X)-p2X]  =0  (B-2) 


where  the  values  X  and  7  contain  0 


X  =  sin2  0 

7  =  0-0 


-\ 


( B— 3) 


For  given  values  of  the  velocity  V  and  of  the  properties 
p  ,  G  ,  v  of  the  material,  Eq.  (B-2)  defines  0  implicitly 
as  function  of  Q  and  (3  ,  which  describe  the  state  of  stress 
at  the  point  considered.  Substitution  of  Eqs.  'B-3)  leads  to 
a  quartic  in  tan  0  ,  so  that  the  direct  investigation  of  the 
character  of  its  roots  is  an  unpleasant  prospect.  An  alter¬ 
native  approach  will  be  used  here. 

Consider  first  the  properties  of  the  roots  X  of  Eq.  (B-2) 
as  functions  of  (3  and  7  .  Equation  (B-2)  is  quadratic  in  X: 
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C  =  4(l-2v) ( 3+P2)x2  -  4[(4-5v)  +  (l-v)p2  +  (l-2v)0  cos  27]X  + 


P2cos2  27  +  2(l-2v){3  cos  27  +  5  -  4v  =  0 


(B-4) 


To  prove  that  this  equation  must  have  two  positive  roots  X, 
it  is  noted  that  examination  of  the  X-independent  term  shows 
C  0  for  X  =0  .  The  coefficient  of  X2  in  the  range 
0(v<  l/2  '  being  positive,  C  is  also  positive  for  suffi¬ 

ciently  large  values  of  X.  The  proof  is  completed  by  the 
demonstration  in  the  next  paragraph  that  there  is  a  positive 
value  X  for  which  c  <(  0 

It  is  proved  in  [6]  in  the  transient  case  that  one  pair 
of  characteristic  velocities  is  always  above,  one  below  the 
speed  cc  of  elastic  shear  waves.  The  equivalent  holds  here 
for  X  and  is  shown  by  substitution  of  the  appropriate  value 

Xs  =  1/2  (B-5 ) 


in  Eq.  (B-2) , 


*) 

The  value  v  =  l/2  need  not  be  considered  as  pointed  out 
in  Appendix  A. 
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c  lx=  1/2  =  P2tcos2  27  "  1]  (B_6) 

Excepting  the  cases  3=0  or  cos  27  =  ±1  when  C  vanishes, 
Eq.  (B-6)  gives  negative  values  of  C,  so  that  there  are  two 
distinct  real  roots  0  <(  X_^<(l/2  ,  X^>l/2  .  The  special 

cases  3=0  and  cos  27  =  ±  1  require  consideration  of  the 
discriminant  D  of  Eq.  (B-4) , 

D  =  (l-2v)32(3+32)  sin2  27  +  [ ( l-2v)3  cos  27  +  1  +  v]2  + 

v32[ v32  +2(l-2v)3  cos  27  +  2(l+v)]  (b-7) 

If  3  =  0  the  value  of  D  becomes  (l  +  v  )2  7^  0  so  that 
X  =  l/2  is  a  root  but  not  a  double  one.  If  cos  27  = ±  1  , 
sin  27  =  0  ,  let 

3=3  cos  27  (B-8) 

and  p2=  32  cos2  27  =  32  .  The  discriminant  (B-7)  is  then 

a  function  of  3  and  v  only,  and  of  the  surprisingly  simple 


form 


COS  27  =  ±  1 


=  [(P-l)2v  +  p  +  1] 


(B— 9) 


The  value  X  =  l/2  will  therefore  be  a  double  root  if 


cos2  27  =  1,  and  v 


B_cos__27_ji_l__ 
(f3  cos  27  -  l)2 


(B-10) 


Allowing  for  the  special  cases,  it  has  been  shown  so  far 
that  for  given  values  {3  and  7  ,  Eq.  (B-2) ,  has  only  real 
positive  roots  X, 

0  <  Xi  <  V2  <  X2  (B-ll) 

The  roots  are  equal,  X1  =  X£  =1/2  if  Eqs.  (B-10)  hold. 

Studying  x(|3,7)  as  a  function  of  7  alone,  relative 
maxima  and  minima  of  X  for  any  given  value  of  (3  must  satisfy 
the  equation 

dC  _ 
cS7 

This  will  be 
(B-12)  gives 


-4p  sin  27 [  (3  cos  27  +  (l-2v)  (1-2X)]  =  0  (B-12) 


used  to  find  an  upper  bound  on  X  .  Equation 

2 

three  possibilities: 


-66- 


{3  cos  2y  =  ~  (l-2v)  (l-2x) ,  or  sin  2y  =  0,  or  6  =  0 
Case  a.  When  £  cos  2y  =  -  (l-2v)(l-2x)  ,  Eq.  (B-2)  becomes 

Ca  =  [ i+(l-2v) (l-2X) ] [ (i+v) (l-2x)  -  p2X]  =  0  (B-13) 

which  gives  either 


X 

a 


i-v 

l-2v 


independent  of  P  ,  or 


(B-14) 


(B-15 ) 


The  maximum  of  the  latter  value,  which  occurs  for  (3  =  0 
is  l/2  .  Xp  being  larger  than  l/2  the  result  (B-15)  never 
controls . 


Case  b.  sin  2y  =  0  .  In  this  case  cos  2y  -  ±1  and  the  value 

of  C  is 


Cb  =  [±£  +(l-2v) (l-2x) ] 2  +  2[l+(l-2v)(l-2X)] 


[ (l+v) (1-2X)  -  P2X] 


(B-16) 


Substituting  the  maximum  value  of  X  found  in  Case  a,  from 


Eq.  (B-14)  into  Eq.  (B-16)  one  finds 


X 


-  1-v 

a"  l-2v 


[±P  "  l]2  >  0 


(B-17 ) 


It  was  previously  demonstrated  that  the  value  of  C  is  not 

positive  for  X  =  l/2  and  that  there  can  be  only  one  root 

X  /  1/2  •  The  value  c,  for  X  =  X  being  positive,  the  only 
'  b  a 

root  of  Eq.  (B-16)  for  which  x^  /  l/2  must  be  less  than 

1—v 

X  =  .  Case  b  cannot  furnish  the  upper  bound, 

a  l-2v  ^ 

Case  c.  In  this  case  p  =  0  and 


Cc  =  (l-2X)[3(l-2v)(l-2x)  +  2(l+v) ]  (B-18) 


This  equation  gives  Xc=  1/2  /  which  is  not  a  maximum,  and 


5-4  v 
6(1-2v) 


(B-19 ) 


Comoarison  with  the  value  of  X  =  ,  Eq.  (B-14) ,  indicates 

a  T? 

that  X  /  Xp  regardless  of  the  value  of  v  . 

While  the  controlling  maximum,  Eq.  (B-14) ,  was  derived 
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by  a  search  for  relative  maxima  of  X  for  given  values  of  ft, 
the  largest  value  of  X  found  being  independent  of  p  ,  it  is 
the  absolute  upper  bound.  For  any  selected  pair  of  values  (3 
and  7  ,  Eq.  (B-2) ,  has  therefore  two  positive  roots,  subject 
to  the  following  bounds: 

0  <  X1  <  1/2  ,  1/2  <  X2  <  Xp  =  (B-20) 

Changing  from  the  values  X  to  the  more  pertinent  values 
0  ,  Eq.  (B-3)  gives 


sin  0 


(B-21) 


The  upper  bound  in  Eq.  (B-20)  ,  X  ,  assures  that  Eq. 

(B-21)  gives  real  values  0  for  superseismic  velocities 


V  > 


/2  ( 1-v )  G 


(B-22 ) 


Restricting  the  range  of  0  to  the  meaningful  one,  0  <(  0  tt 
it  has  therefore  been  demonstrated  that  Eq.  (B-2)  has  two  pairs 
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of  real  roots  0^  }  7r  -  0^  ,  and  no  others,  for  any  pair  of 

values  7  and  (3  .  Double  roots  0  =  0  occur  when  Eqs .  (B-10) 

S 

apply.  The  roots  are  bounded 


0  <  0X  '  7T  -  0  , 

"■ O 


0p  K.  0  <  0 

P  —  3  — 


7T-0  <(  0  <(77-0 

S  —  2  _ 


P 


0  <  0  <  7T 

S  —  4 


7 


}  (B-23) 


The  final  step  may  now  be  taken,  where  the  roots  0  of 
Eq.  (B-2)  are  considered  as  functions  of  (3  and  6  .  For  a 
given  value  of  (3  ,  any  root  (b.  satisfying  Eq.  (B-2)  for  a 
specified  value  7  ,  gives  a  combination  0  and  6  =  0  -  7 
which  satisfies  Eq.  (B-2) .  It  will  now  be  demonstrated  that 
variation  of  7  leads  to  solutions  for  any  value  of  6  . 

Figure  B-l  shows  a  typical  plot  of  0(7)  for  a  fixed 
value  of  (3  .  There  are  four  curves,  bounded  by  horizontal 
lines  representing  the  bounds,  Eq.  (B-23).  The  value  7  appears 
in  Eq.  (B-21)  only  in  the  form  cos  2y  ,  so  that  0(7)  must 
be  periodic,  0(7)  =  0(7 ±  7r)  .  Each  of  the  four  separate 
branches  must  be  continuous  as  function  of  7  .  (There  can  be 
no  discontinuity  because  the  standard  form  of  solution  of  the 
quadratic  equation  (B-4)  for  X  cannot  lead  to  a  discontinuity 
if  X  is  always  real.)  The  values  0  which  pertain  to  a  pre¬ 
scribed  value  of  9  appear  in  Fig.  B-l  as  intersection  points 
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of  a  straight  line,  0  =  7  +  ©  ,  with  the  curves  0  =  0(7)  . 
There  is  obviously  at  least  one  value  0  in  each  interval 
between  the  five  horizontal  bounding  lines,  for  a  total  of 
four,  except  for  the  possibility  of  double  roots  which  may 
occur  for  certain  values  of  9  if  the  second  Eq.  (B-10)  is 
satisfied.  The  purely  qualitative  Fig.  B-l  does  not  prove 
that  the  straight  line  0  =  7  +  0  cannot  intersect  the  same 
branch  more  than  once.  This  is,  however,  impossible  because 
Eq.  (B-2)  is  a  quartic  having  at  most  four  roots.  The  fact 
that  Eq.  (B-2)  as  function  of  the  state  of  stress,  9  and  (3 
is  a  quartic,  not  a  bi-quadratic  indicates  also  that  the  char¬ 
acteristic  directions  are  not  symmetric,  i.e.  in  general 

0  7*-  7T  — 0  J0  7T  —  0 

1  4  2  3 

It  has  therefore  been  demonstrated  that,  including  double 
roots,  there  are  always  four  characteristic  directions  7r)>  0^^>  0 
for  any  state  of  stress,  the  values  0^  being  subject  to  the 
inequalities  (B-23) .  Due  to  the  asymmetry  mentioned  in  the 
last  paragraph,  only  one  double  root  0  )>  0  may  occur  for  any 
state  of  stress,  the  other  two  roots  remain  different.  In 
addition  there  is  a  triple  root  0=0 

In  final  conclusion  it  has  been  proved  that  all  seven 
characteristic  directions  are  real,  and  the  system  is  there¬ 
fore  hyperbolic.  The  occurrence  of  multiple  roots  does  not 
create  any  difficulty  affecting  the  nature  of  the  differential 
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equations  [7]*-  While  the  details  are  not  given  here,  it  was 
found  that  along  a  double  (triple)  characteristic  two  (three) 
independent  compatibility  relations  apply. 
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Result  for  v  =  0,25,  V  *  1.25  cp  ,  pQ  =  3.58k 
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Fig.  14c 
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